Prom self-consistent covariant effective field theories to their Galilean-invariant 

counterparts 
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We discuss how to obtain the nonrelativistic hmit of a self-consistent relativistic effective field 
theory for dynamic problems. It is shown that the standard v/c expansions yields Galilean invariance 
only to first order in v/c, whereas second order is required to obtain important contributions such 
as the spin-orbit force. We propose a modified procedure which is a mapping rather than a strict 
v/c expansion. 
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I. INTRODUCTION 

Since the development of the special theory of relativ- 
ity classical physics is viewed as its nonrelativistic limit 
for C; the velocity of light, going to infinity (or corre- 
spondingly, considering all relevant velocities to be much 
smaller than c). In Quantum Mechanics, Dirac's theory 
of the electron carrying spin 1/2 has as its nonrelativistic 
limit the Pauli equation. 

Since then, the problem in obtaining a nonrelativis- 
tic limit of a theory has enjoyed permanent interest. 
Many physical systems are on the borderline between 
relativistic and classical physics, determined by the ve- 
locity and/or the size of the system. Prominent exam- 
ples are heavy atoms and nuclei, where the spin-orbit 
force is a relativistic effect and shows that, although a 
nonrelativistic description is in general easier to handle, 
relativistic effects cannot always be neglected. Under- 
standing the transition from relativistic to nonrelativistic 
physics while maintaining the traces of relativistic effects, 
is therefore of utmost interest. 

One of the first, and still most prominent, attempts to 
handle this transition in a systematic way was the Foldy- 
Wouthysen transformation [l| using a canonical trans- 
formation to obtain two equations with two components, 
one of which becomes the Pauli equation in the nonrela- 
tivistic limit. An alternative method is related to group 
contractions [1], where one discusses under which condi- 
tions a group can be contracted to another one, involving 
non singular transformations. In particular, the Lorentz 
group 50(3, 1) can be contracted to the Galilei group, 
taking the limit c —s- oo and neglecting corrections of 
the order of (f )^. On the level of group generators this 
contraction procedure is direct, but on that of represen- 
tations it is not as trivial i. As we shall see, this is 
the case in self-consistent effective field theories involv- 
ing fields with a spin i representation. In recent years, 
several groups discussed similar problems ^, i^, i^] for the 
Maxwell and Dirac equations. They showed that several 
nonrelativistic limits may exist in these cases. 



In the present paper, we focus on the nonrelativistic 
limit of covariant nonlinear self-consistent effective field 
theories, i.e. those based on density functionals. These 
cases are to be distinguished from other effective field 
theories which rely on a straighforward expansion (see 
e.g. @) and for which the following consideration do 
not apply. Realizations of self-consistent field theories in- 
clude e.g. the model by Duerr Q, Heisenberg's nonlinear 
spinor theory Q, the Nambu-Jona-Lasinio (NJL) model 
0, and effective c/)^ theories [l3|. In nuclear physics, 
the relativistic mean-field (RMF) model ll|, 113 ^^i^ the 
Skyrme-Hartree-Fock (SHF) approach [l2| are typical ex- 
amples. 

In this manuscript, we will discuss the link be- 
tween such a non-linear relativistic theory and its non- 
relativistic counterpart. It is found that a straightfor- 
ward nonrelativistic reduction of a covariant ansatz up to 
(v/c)^ yields a result which violates Galileian invariance. 
We will develop and justify a non-relativistic mapping 
going up to {v/cY that leads to the correct Galilean- 
invariant counterpart. 



II. NONRELATIVISTIC REDUCTION 

A. The goal 

A relativistic effective field theory expresses the con- 
figuration in terms of Dirac four-spinor wavefunctions V'q 
for each state a. For the further developments, it is use- 
ful to express it explicitly through upper and lower two- 
spinor components as 




(1) 



For means of simplicity and lucidity we choose a simple 
covariant self-consistent effective field theory involving 
point couplings between its degrees of freedom, i.e., the 
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4-componcnt spinors. It reads 



with 



a a 

Q a. 



(2) 



(3a) 
(3b) 
(3c) 



with 7^ = (70 , 7) the four-veetor of Dirae matriees [ij] • 
The relativistic functional appears simple because the ki- 
netic and spin-orbit terms are implicit in the scalar and 
vector densities, as we shall see. For simplicity of no- 
tation, we will drop the index label a in the following, 
identifying, e.g., the scalar density with Qs = ip^l' and 
similarly for all other densities and currents. 

A key point in such approaches is, of course, correct 
normalization of the densities and consequently the wave- 
functions. This implies 



d'^r-!/i7o'!/' 



(4) 



to guarantee invariance under Lorentz transformations. 
This may be surprising, since it involves the zeroth com- 
ponent of a four vector, but is inevitable to counter the 
relativistic contraction of the purely spatial volume ele- 
ment d'^r [l^l- We will see that this is the key problem 
with straightforward expansions and at the same time 
the key to the solution. 

The goal is now to obtain a nonrelativistic functional 
based on the densities © from the relativistic parent 
functional ([2]) in a nonrelativistic limit. Solutions in the 
positive-energy branch (particle-like) arc distinguished 
by a dominance of (p^^^ over (p'^'^\ The strategy is thus to 
eliminate ip^^^ and identify the upper component with the 



classical two-spinor wavcfunction, ip^^^ ^ 
should obey the normalization condition 



dV|^(^')| =1. 



ip^'^^\ which 



(5) 



Wc can thus express ip in terms of i^''^''. Inserting that 
into the coupling Lagrangian density ^ should produce 
the desired limit. 

In the course of obtaining the nonrelativistic limit we 
expect to formulate the nonrelativistic counterpart in 



terms of the following densities and currents: 



(6a) 



-^E['^'"'^(^x^)^''- (^X'^^-V'P^.'], (6b) 

a 

-^E[^'-'^^V'^"'-(^^«)Va], (6C) 

(6d) 



where p, t and J are time-even while j and cr arc time- 
odd. Note that time-even and time-odd terms appear 
in particular combinations, a feature which is crucial to 
render the functional invariant under Galilean transfor- 
mations llJl. 



B. Problems with ii/c expansions 

Nonrelativistic limits are usually obtained by straight- 
forward expansion in orders v/c or p/m, respectively, 
e.g. in the Foldy-Wouthu yse n transformation ([3], for 
the nuclear case see [ll|. 111]). We briefly review the 
steps from [llj. One starts from the relativistic equa- 
tions of motion, for the present model derived from the 
Lagrangian ^ 

= (i7^a^-m + 5 + 7^T/^)Va, (7a) 

5 = -CsQs, Vf_, = C„£)p. (7b) 

with the self-consistent scalar and vector potentials S and 
V^, the latter decomposing as 

= (l/o,-V) = (c„go,-c„e). (8) 

We insert the decomposition ^ and solve the lower- 
component equation for (^'•'^■' . Keeping only terms up 
to order p/m yields 



Bo = 



Bo^-(p-V)^("), 
1 _ 1 

2m + 5 - Vo ~ 2m' 



(9a) 
(9b) 

1/m. ig- 



where p = — iV. The approximation Bq 
nores the density dependence in Bq. It suffices for the 
present studies. The form ([9]) violates the normalization 
(III) to second order in p/m. The procedure of [lH restores 
(ortho-)normalization at operator level by first imposing 
the relativistic normalization ^ up to (p/m)'^ 



drrp 



(u)t 



r = i?2(^.(p_v)) 



i + r 
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■•o y" ' yy ~ ) ) ■ (10) 

We introduce ip^'^^'' so as to recover the nonrelativistic 
normalization ([5]). Thus we identify 
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1/2 


i + t 
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We expand in second order of p/m and obtain 



2 . 



(cl) 



(11) 



which, together with relation ([9a|) . provides a complete 
description of ^/^ in terms of i^s^'^'-' up to order p/m. Using 
the relation (o--A)(<t-B) = A-B — iA- (it x B) yields 

^ = (p-V)2 + i(p-V).(^x(p-V)). (12) 

We insert eqs. ©, pT|) and into expressions ([3]) for 
the relativistic densities and finally obtain up to order 

{p/mf 



Qs = p-'^Bl T-V-J-c„£.- (2j + Vxcr) 



go 



2-Bo(j -CyQp + -^Vxcr 



(13a) 
(13b) 



The same relations are obtained when going through 
the Foldy-Wouthuysen transformation up to second or- 
der. Note that the spatial part of the vector density q 
is at least of first order such that the correction from 
T would be of third order and is discarded. The result 
is correct for stationary states where all time-odd den- 
sities and fields vanish, i.e. where j = 0, <t = 0, and 
V = 0. The expressions (|13p when inserted into the 
Lagrangian ^ produce a serious defect: the emerging 
nonrelativistic Lagrangian is not Galilei invariant. The 
Lorentz invariant scalar density ps does not translate to 
a Galileian invariant expression and the same happens 
for the combination — ■ Correct expressions should 
form the Galileian invariant combinations pr — and 
pV • J + j • (V X cr) as we will see later. 

In order to elucidate the problem, we consider the 
transformation properties for the simple case of an ex- 
plicit boost of the whole system. Let us start with a 
well-checked situation, a stationary state for which j = 
and <T = 0. We boost the system by a constant velocity 
u (in units of c). All quantities in the boosted system 
will be distinguished by a prime. The normalization (|4]) 
becomes in the boosted frame 



1 = j dV £-0 = y" dVVl - u2 



£'0 



(14) 



Note the transformation of the volume element exactly 
countering that of the density qq. The volume dilatation 
factor is negligible to order but contributes in second 
order. The above nonrelativistic expansion to second or- 
der had violating terms at that order. The example shows 
that the mistake lies in neglecting a second order correc- 
tion of the volume element in the normalization condi- 
tion. 



C. Map instead of expansion 

The previous discussion shows that a straightforward 
nonrelativistic expansion with all kinetic contributions is 



consistent only up to first order p/m (or boost velocity u, 
respectively), while the crucial relativistic corrections to 
a classical Schrodinger equation emerge from second or- 
der terms, namely spin-orbit coupling and effective-mass 
terms. These require a special handling of the normal- 
ization condition like in the example of Eq. (fT4|l . We 
thus leave the straightforward paths oi p/m expansion 
to now aim at a generalized mapping of the relativistic 
functional into a nonrelativistic one, trying to incorpo- 
rate all second-order effects. 

The key point is accounting for the relativistic volume- 
element compression as in Eq. (fT4|) . To deal with arbi- 
trary nonstationary situations we need to allow a boost 
velocity field. Thus we modify the normalization condi- 
tion ([4]) to display explicitly the compression factor with 
respect to the local boost velocity, which in turn is ex- 
pressed in terms of the classical densities and currents: 



1 = j d'Vv'l-u2(r) 



(15a) 



u(r) = -^^^(j-Vp+ivxa). (15b) 
qq p \ 2 / 

The expansion (fTTjl thus is slightly modified to the map 



1 - -r 
2 . 



^(^')(l + ^u^) 



1^ 1 9 

1 - -r + -u^ 
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(cl) 



(15c) 



Things now proceed as in section lll Bl but the term cx 
cancels the unwanted one in and adds a desired one 
in go- This now leads to the consistent result 



pr-j2 -pV-J-j-(Vxcr) 



-(Vxo-)^ 



1. 



Qq = P + 



Q = 2Bo(j- Vp+-Vxcr 



(16a) 
(16b) 
(16c) 



The scalar density shows the wanted Galilean-invariant 
combinations and the vector density reproduces the cor- 
rect invariance property, q^q'^ ^ Pq ~ ~ p'^ ^ up to 
terms of second order, of course. We thus insert the 
mapped scalar and vector densities p6p into the interac- 
tion Lagrangian density ^ , getting (up to second order) 



PT-3 



-(pV-J+j-(Vx^))--(Vx^)2 



(17) 



That result is manifestly Galilean invariant. 

It has the same form as the basic version of the Skyrme 
Hamiltonian density [l^ that is being employed for the 
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description of finite nuelei. This comes as no surprise, 
since the covariant Lagrangian that we started with, Eq. 

consists of the basic terms of the Lagrangian of the 
point-coupUng variant of the RMF model for nuclear 
structure, RMF-PC [H]. Both SHF and RMF-PC are 
formulated in terms of point couplings of spinors and thus 
display this close relationship. This relation between co- 
variant models and their Galilean-invariant counterparts 
is of importance when one analyzes, for example, spin 
excitation mechanisms in nonrelativistic time-dependent 
Hartree Fock employing the Skyrme functional • 

Compared with the Skyrme functional, however, there 
is one additional term oc (Vxcr)^ which is an allowed term 
in the Skyrme functional, but usually neglected. There is, 
however, no gradient term oc pAp which is mandatory for 
the description of finite-size systems. That is no surprise 
because we had started from a simplified Lagrangian 
without gradient terms. The more complete model would 
also include terms such as PsAp^ and p^Ap^, whose ex- 
pansion proceeds quite similarly and in the nonrelativis- 
tic limit yields gradient terms. There is a subtle dif- 
ference, though: the nonrelativistic mapping would also 
produce gradient kinetic terms like pAr, which are ne- 
glected assuming that the gradient correction as such is 
small and second order relativistic corrections to it arc 
negligible. The counter argument is that there are two 
quite different notions of smallness involved here that 
may not be combined. 



as to the current and spin densities in the nonrelativis- 
tic domain. Previous derivations rarely studied the full 
dynamical case. 

The procedure started out with a straightforward v/c 
expansion, which encountered inconsistencies, because 
the spin-orbit term appears only in second order of v/c 
while Galilean invariance comes out correctly only to first 
order. The key finding is that a strictly nonrelativistic 
theory is not easily compatible with the appearance of a 
spin-orbit term. To be more precise: for relativistic effec- 
tive field theories, based on density functionals, it is not 
possible to derive a sufficiently complete nonrelativistic 
theory by mere expansion and order counting. Instead 
the more general concept of a nonrelativistic mapping 
is needed, namely to derive an effective nonrelativistic 
theory which includes as many features of the given rela- 
tivistic theory as desired. Starting from the simple con- 
sideration of Lorentz contraction of the spatial volume 
element, we have derived such a mapping for a covariant 
self-consistent model. This mapping manages to provide 
a manifestly Galilean invariant theory which correctly in- 
corporates the spin-orbit and effective-mass terms and, if 
one starts with the RMF-PC model in nuclear physics, 
merges into the widely used Skyrme-Hartree-Fock ap- 
proach when neglecting the involved density dependences 
of the spin-orbit and effective-mass term. Extensions of 
the scheme developed here are in progress. 



III. CONCLUSION 

We have studied the nonrelativistic limit of a self- 
consistent relativistic theory with the aim of recovering 
a minimum of relativistic effects, the spin-orbit force, to- 
gether with Galilean invariance in the resulting nonrel- 
ativistic theory. Note that Galilean invariance requires 
keeping all time-odd terms, which play a crucial role in 
the formulation of dynamics. This applies to the spa- 
tial components of the relativistic vector density as well 
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